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I. INTRODUCTION 



A. The top quark 



The standard model (SM) of particle physics describes 
the elementary particles and their interactions. The top 
quark (t) has a special place in the hierarchy of particles 
because it is far more massive than any of the other 
fundamental objects. It is the up-type quark, partnered 
with the down-type bottom quark (b), forming the third 
generation of quarks which was predicted by Kobayashi 
and Maskawa in 1973 [lfto accomodate CP violation 
in neutral kaon decays Q- At particle colliders the 
top quark is produced mainly in top-antitop (tt) pairs. 
First evidence of top quark production was reported by 
the CDF collaboration Q and the top quark was first 
observed in this production mode by the CDF [|| and 
DO [j| collaborations at the Tevatron proton-antiproton 
collider. Since then, great efforts have been focused 
on measuring its properties with ever higher precision. 
In addition to its large mass (to*), the top quark is 
also singular because it decays before it can hadronizc: 
there are no mesons or baryons containing valence top 
quarks. The top quark decays almost exclusively to a 
W boson and a b quark, with the fraction determined 
by the near-unity value of the CKM quark mixing 
matrix [HH element V tb (« 0.9992) Q. Its other decays 
are limited by the small values of V ts « 0.0387 and 
Vtd ~ 0.0084 [31, assuming three-family unitarity of the 
CKM matrix. The W boson decays to a charged lepton 
and its associated neutrino, or to a quark- antiquark pair, 
and the final states of tt events are thus characterized as: 
"lepton+jets'^ (tt -> £+vbqq'b and qq'b£-Db); "alljets" 
(tt ->■ qcfbqq'b), and "dileptons" (tt ->■ vbt 9b) . In 
this notation the charged lepton £ represents an electron 
or muon, and q is a first- or second-generation quark. The 
W boson also decays to a r lepton and a r neutrino. If the 
r decays to an electron or muon, the event contributes to 
the lepton categories, and if the r decays into hadrons, 
it contributes to the lepton+jets or alljets categories. A 
fourth category labelled "^r+jets" is used to measure mt 
when there are jets and a large imbalance in transverse 
momentum in the event ($t), but no identified lepton. 
It comprises tt — > T + vbr~ Db, T + ubqq'b, and qq'bT~vb 
final states, accounting for 40% of the tt signal events in 
the ^T+jcts category, or l + vbqq'b, qq'b£~Db, where the 
electron or muon are not reconstructed, accounting for 
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60% of the tt signal in this category. Additional contri- 
butions to arise from the neutrino(s) produced in r 
decays. 

In dilepton events, there are typically two jets from 
the two b quarks, one from each top quark decay. In 
lepton+jcts events, there are typically four jets, including 
two b jets and two light-quark jets from VF-boson decay. 
Alljets events most often contain six jets, the two b jets 
and four light-quark jets. The $r+jets events usually 
have four or five jets. Additional gluon or quark jets can 
arise due to radiation from initial or final state colored 
particles, including the top quarks. About 23% of the 
tt events have an extra jet with sufficient energy to pass 
the selection criteria and about 5% of the events have 
two additional jets. These extra jets complicate the 
measurement of mt and degrade its resolution. Figure Q] 
illustrates leading-order (LO) production of tt events at 
the Fcrmilab Tevatron Collider, and Fig. [2] shows the 
relevant tt decay modes. 
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FIG. 1: Examples of tree Feynman diagrams for tt 
production. At the Tevatron collider, the qq channel 
contributes 81% to the total tt inclusive cross section and 
the gg channel the remaining 19% 




e, \i~ t, q 



FIG. 2: Leading-order Feynman diagram for tt decay. The 
dilepton modes (ee, e/i, /i/x) have a combined branching 
fraction of ~ 4%, the electron+jets and muon+jets modes 
combined correspond to ~ 30%, and the alljets mode has a 
branching fraction of « 46%. The r modes are shared among 
the ^?T+jets and the other channels in the analyses. 



B. Top-quark mass origin and definitions 

One of the fundamental properties of an elementary 
particle is its mass. In the SM, fermions acquire mass 
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through interactions with the Higgs field 0. Absolute 
values of these masses are not predicted by the SM. In 
theoretical calculations, a particle's mass can be defined 
in more than one way, and it depends on how higher- 
order terms in perturbative quantum chromodynamics 
(QCD) calculations are renormalized. In the modified 
minimal subtraction scheme (MS), for example, the 
mass definition reflects short-distance effects, whereas in 
the pole-mass scheme the mass definition reflects long- 
distance effects [ll[ . The concept of the pole mass is not 
well defined since color confinement does not provide S- 
matrix poles at m = m t [l2j . Direct mass measurements 
that are inputs to the combination described in this 
paper rely on Monte Carlo (MC) generators to extract 
mt- Hence the measured mass corresponds in fact to the 
mass parameter in the MC. Work is proceeding to address 
the exact difference between the measured mass and the 
pole mass, as presented for example in Appendix C of 
Ref. |13j . One alternative way to address this problem 
is to extract mt from a measurement of the tt cross 
section [l4|- The DO collaboration has recently shown 
that the directly measured mass of the top quark is closer 
to the pole mass extracted from a measurement of the tt 
cross section than to an MS mass extracted in a similar 
way [l4[. Hence, within the precision of theory and data, 
the directly measured mt is best interpreted as the top- 
quark pole mass. 

CPT invariance predicts that a particle and its 
antiparticle partner have the same mass. This has been 
checked for the top quark by the DO, CDF, and CMS 
collaborations and the masses are found to hold within 
the measurement uncertainties, with Am t = m± — fflf = 
0.84 ± 1.87 GeV Am, = -3.3 ± 1.7 GeV [H], and 
Am t = -0.44±0.53 GeV [3], respectively. Thus, the top 
quark mass combination in this paper assumes mt = mi. 

C. Predictions based on the top-quark mass 

The internal consistency of the SM can be tested by 
using different observables to predict the values of others, 
and then compare the expectations with their measured 
values. For example, the relation between the mass 
of the W boson {My/) and sin 2 6\y (the electroweak 
mixing angle) includes higher-order radiative corrections 
involving m t , hence the smaller the uncertainty on the 
measured m t , the stronger is the test of consistency. 

Since 1997, the LEP Electroweak Working Group 
(LEPEWWG) has used the the observed top-quark and 
the W boson masses and other precision electroweak 
variables to extract constrains on the Higgs boson mass 
(M H ) in the SM QJ]. This has been extended to the 
minimal supersymmetric standard model (MSSM) [l9[ 
and the GFITTER collaboration has applied the 
technique to set limits on a wide variety of theories 
beyond the SM |20| . Figure [3^, shows the combined 
constraint due to M\y and m t (as of March 2012) on the 
Higgs boson mass. Figure ISb shows the constraint from 
Mw and m* separately (as of March 2012) on the Higgs 



boson mass, and a global constraint originating from all 
the other electroweak variables, showing the importance 
of the Mw and m t variables to contrain the Higgs boson 
mass. 
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FIG. 3: Panel (a) shows constraints from LEP and Tevatron 
measurements of Mw and m t (Tevatron only) on Mh within 
the SM. The regions in mass of the Higgs boson still allowed 
after the direct searches at LEP, Tevatron and LHC are also 
shown. Panel (b), from Ref. [2(j, indicates with the large 
countours (blue) the contraints on the Higgs boson, from 
global fits to electroweak data without including the direct 
measurements of Mw and mt from the Tevatron. 

D. History of measurement of m t 

Before 1995, global fits to electroweak data from the 
CERN and SLAC e+e" colliders (LEP and SLC) and 
from other experiments produced estimates of mt that 
ranged from w 90 GeV to « 190 GeV 0. At the 
time of the first observation of the top quark in 1995, 
the fits indicated a mass close to the current Tevatron 
value of mt, but with an uncertainty of « ±10% and an 
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assumption of 300 GcV mass of the Higgs boson [231 . 
CDF measured m t = 176 ± 8(stat) ± 10(syst) GeV~|l 
(total uncertainty of 7%) and DO measured m t = 
199±2?(stat) ± 22(syst) GeV [H] (total uncertainty of 
15%). 

Since then, the CDF and DO collaborations have 
developed many novel measurement techniques, 
and published nearly 50 journal papers on their 
measurements of mt. Recently, the CMS collabo- 
ration at the Large Hadron Collider (LHC) published 
a measurement using 102 dilepton events [23[ and 
finds m t = 175.5 ± 4.6(stat) ± 4.6(syst) GeV 
(total uncertainty of 3.7%). The ATLAS collabo- 
ration at the LHC has submitted a measurement of 
ra t = 174.5 ± 0.6 ± 2.3 GeV (total uncertainty of 1.4%) 
using nearly 12,000 lepton+jets events (24|. The most 
precise measurements from the Tevatron in a single 
decay channel use lepton+jets events and a matrix 
element method as introduced in Ref. [25| and an in-situ 
calibration of the jet energy scale. CDF's matrix element 
measurement [26| uses 5.6 fb _1 of integrated luminosity 
to find m t = 173.00 ± 0.65(stat) ± 1.06(syst) GeV 
(total uncertainty of 0.72%). DO's measurement [27j 
uses 3.6 fb _1 of integrated luminosity to obtain 
m t = 174.94 ± 0.83(stat) ± 1.24(syst) GeV (total 
uncertainty of 0.85%). Figure 0] shows the publication 
history of the direct measurements of mt at the Tevatron. 

E. Overview of mass measurements 

This paper reports on the combination of previously 
published measurements of m t . Details of the analyses 
are therefore not repeated as this information is available 
in recent reviews [28|], as well as in the publications of 
each of the results. We will, however, summarize the 
basic techniques used for the measurements. 

The cross section for ti production in proton- 
antiproton (pp) interactions at 1.96 TcV is 7.2 pb (29l . 
[30I ] - The mean transverse momentum (pt) of the ti 
system at parton level is k, 20 GcV, which is attributed 
to initial-state radiation (i.e., gluon emission). The 
mean transverse momentum of the top quarks at parton 
level is « 95 GeV 13111 . Top quarks have a lifetime of 
« 0.3 x 10~ 24 s (32l. |33| . which is an order of magnitude 
smaller than the time scale for parton evolution and 
hadronization. Hence, when top quarks decay, they 
transfer their kinematic characteristics to the W boson 
and b quark, and the measured energy-momentum four- 
vectors of the final state particles can be used to 
reconstruct the mass of the top quark, except for the 
presence of initial or final-state radiation. 

In alljets events, the four- vector of every jet emerging 
from quarks can be reconstructed, but neutrinos emitted 
in semileptonic decays of b quarks and jet energy 
resolution effects will lead to lost energy. In lepton+jets 
events, the momentum of the neutrino from the W — > 
£f£ decay is not detected. The transverse component 
can be inferred from the negative of the vector sum 
of all transverse momenta of particles detected in the 



calorimeter and muon detectors. We estimate the 
longitudinal momentum of i>i by constraining the mass 
of the charged lepton and neutrino system to the world 
average value of M\y 0j- We also use My/ to choose 
the two light jets from W —> qq' decay, and use that 
information for an in-situ calibration of jet energies. In 
dilepton events, the analysis is more complicated because 
there are two final-state neutrinos from the leptonic 
decays of both W bosons. Therefore, the longitudinal 
and transverse momentum components of the neutrinos 
can not be determined without the application of more 
sophisticated tools. These involve assuming a value for 
m t to solve the event kinematics and assigning a weight 
to each mt hypothesis to determine the most likely value 
of m t consistent with the hypothesis that the event is a 
ti event. 

A major issue in ti final state reconstruction is 
the correct mapping of the reconstructed objects to 
the partons from the decays of the top quark and 
W boson. The problem arises because often the 
jet charge and flavor cannot be uniquely determined. 
This creates combinatorial ambiguities in the ti event 
reconstruction which vary from 90 possible jet-to-parton 
assignments for the alljets final state to 2 in the 
dilepton channel. In the lepton+jets and dilepton final 
states, additional ambiguities may arise from multiple 
kinematical solutions for the longitudinal component of 
the neutrino momentum. 

Two methods are used to measure the value of mt- 
In the first method, the reconstructed mass distribution 
in data, or a variable correlated with m t , such as 
the decay length of the B hadron or the transverse 
momentum of a lepton, is compared to template distri- 
butions composed of contributions from background 
and simulation of ti events. One template is used to 
represent background and another for each putative value 
of m t . The second method uses event probabilities 
based on the LO matrix element for the production of 
ti. For each event, a probability is calculated as a 
function of m t that this event is from ti production, as 
based on the corresponding production and decay matrix 
element. Detector resolution is taken into account in 
the calculation of these probabilities through transfer 
functions that correlate parton-levcl energies and their 
measured values. The value of mt is then extracted from 
the joint probability calculated for all selected events, 
based on the probability for signal and background 
(also defined through its matrix element). This method 
produces the most accurate results, but the computations 
are time-consuming. 

F. Combination overview 

This paper describes the combination of statistically- 
independent top-quark mass measurements from the 
Fermilab Tevatron Collider. Measurements are 
independent if they are based on different data sets, 
e.g., from CDF and from DO, or from Tevatron Run I 
(1992-1996) and Run II (2001-2011). They are also 
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independent within one data set if the event selections 
are designed to be exclusive, i.e., no event can pass more 
than one category of selections. At times, more than 
one measurement is published using the same data and 
decay channel. In this situation, the result with smallest 
overall uncertainty is chosen for the combination. Twelve 
measurements are used in the combination described 
here, eight from the CDF collaboration and four from 
DO. These comprise five lepton+jets measurements 
(CDF and DO, Run II and Run I, and a CDF Run II 
result based on the decay length of B hadrons); two 
alljets measurements (CDF Run II and Run I); four 
dilepton measurements (CDF and DO, Run II and 
Run I); and a ^r+jets measurement (CDF Run II). We 
combine these measurements using an analytic method 
called the best linear unbiased estimator (BLUE) [3^ - 
[36| . This technique forms a linear combination of the 
separate unbiased mass measurements to produce the 
best estimate of mt with the smallest uncertainty. This 
procedure follows a series of 11 such mass combinations 
presented in [37l - l47| , updated each year since 2004 as new 
measurements of m t became available. The combination 
presented here is the first to be published in a peer- 
reviewed journal. 

II. INPUTS TO THE COMBINATION 
A. The independent mass measurements 

The mass measurements included in the combination 
are shown in Table [fl [MHIiHIz|. These 12 channels 
are chosen because they are statistically independent, 
which maximizes the improvement in the combination, 
and because enough information is available to separate 



out the components of systematic uncertainty for proper 
treatment in the combination. 

The DO measurement from 2005 in the alljets 
channel (Run I) HH of m, = 178.5 ± 13.7(stat) ± 
7.7(syst) GeV (total uncertainty of 8.8%) is not included 
in the combination because some subcomponents of the 
systematic uncertainty are not available. 

The CDF measurement from Run II based on decay- 
length analysis [57[ differs from the others in that it 
uses the mean decay length of B hadrons in 6-tagged 
lepton+jets events as the mt-sensitivc variable. It is 
independent of energy information in the calorimeter, 
and its main source of systematic uncertainty is 
uncorrelated with the dominant ones from the jet 
energy scale calibration in other measurements. This 
measurement of m t is essentially uncorrelated with 
the higher precision CDF result from the lepton+jets 
channel. The overlap between the data samples used for 
the decay-length method and the lepton+jets sample has 
therefore no effect. 

B. Data 

The data were collected with the CDF f5(| and 
DO [13, HH detectors at the Tevatron pp collider at 
Fermilab between 1992 and 2009. The Tevatron "center- 
of-mass" energy was 1.8 TeV in Run I from 1992 to 
1996 and 1.96 TeV in Run II from 2001. A silicon 
microstrip tracker around the beampipe at the center of 
each detector was used to reconstruct charged-particle 
tracks (only in Run II at DO). Tracks spatially matched 
to calorimeter jets are checked for originating from a 
secondary vertex, or for evidence that they originate 
from decays of long-lived heavy-flavor hadrons containing 
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TABLE I: Top quark mass measurements used as input to determine the combined value of mt from the Tevatron and the 
combined result. 
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b quarks from the decay of top quarks [59l . l62j . Electrons 
and jets produce particle showers in the calorimeters, 
and the collected information is used to measure their 
energies. Muons traverse the calorimeters and outer 
muon detectors that are used to reconstruct their tracks. 
Both CDF and DO have central axial magnetic fields 
in the tracking region (DO only in Run II), in which 
the momenta of charged particles are determined from 
the curvature of their tracks. The CDF magnet has a 
diameter of 3 m and extends 4.8 m along the beamlinc, 
with a field strength of 1.4 T, and the DO magnet 
has a diameter of 1.0 m and length of 2.7 m to fit 
inside the Run I calorimeter with a field strength of 
2.0 T. The CDF detector's larger tracking volume with a 
higher density of measurements gives better transverse- 
momentum resolution for charged-particlc tracks. The 
transverse-momentum resolution is ~ 3.5% at CDF and 
ps 10 % at DO for a muon with p T = 50 GeV. The 
trigger and event-selection criteria depend on the tt 
final states, details of which appear in the publications 
listed in TablclU The experiments collected O(10 14 ) hard 
collisions, from which 7420 events are selected because 
they have the characteristics expected for tt pairs, of 
which ss 56% are expected to be true tt events. 

C. Models for tt signal 

The tt signal in Run I was simulated u sing the 
LO generator herwig [IH with the MRSD [(35f and 
CTEQ4M [H parton distribution functions (PDF) used 
by CDF and DO, respectively. The HERWIG generator 
implements the hard scattering processes qq^ti and 
gg^-tt, adding initial-state and final-state radiation 
through leading- log QCD evolution (6?J ■ The top 
quark and W boson in herwig decay according to the 
branching fractions listed by the Particle Data Group Q , 
and the final-state partons are subsequently fragmented 
into jets. The MC events are then processed through a 
fast simulation or a GEANT model [68| of the detectors 
and then through event reconstruction programs. 



For the tt signal in Run II, CDF uses pythia (6§| 
with the CTEQ5L [z3| PDF, and DO uses the leading- 
log generator ALPGEN [zH with the CTEQ6L1 PDF 
and pythia for parton showering. ALPGEN contains 
more tree-level graphs in higher-order a s than PYTHIA. 
ALPGEN has parton-jet matching (72|, which avoids 
double-counting of partons in overlapping regions of jet 
kinematics. CDF sets the event generation factorization 
and renormalization scales Q 2 to m 2 + p\ + (P 2 + 
PaJ)/2, where p± is the transverse momentum charac- 
terizing the scattering process, and P 2 and P| are the 
virtualities of the incoming partons. DO sets the scales 
to m 2 + (p 2 -), where (p 2 -) is the average of the square of 
transverse momentum of all other light partons produced 
in association with the tt pair. The pythia model 
treats each step of the tt decay chain it— >Wb, W^-lv or 
qq') separately and does not preserve spin correlations. 
ALPGEN uses exact matrix elements for each step and 
thereby correctly describes the spin information of the 
final-state partons. The fragments of the proton and 
antiproton or "underlying event" are added separately to 
each hard collision. CDF uses the "Tune A" settings [73| 
in pythia while DO uses a modified version of the tune. 
Both collaborations use angular ordering for modeling 
parton-showering in PYTHIA, and not pr-ordered models. 
The underlying event is therefore not interleaved with the 
parton showers as in models of color reconnection [74| . 

D. Background models 

In the lepton+jets channel, the dominant background 
is from W^+jets production. Smaller contributions arise 
from multijet events, Z+jets, single top-quark (tqb and 
tb), and diboson production (WW, WZ, and ZZ). The 
alljets channel has mainly multijet events as background. 
The largest background in the dilcpton channel is from 
•Z+jets events, which include Drell-Yan production. 
Backgrounds from diboson production and from events 
with jets identified as leptons are very small in the 
dilcpton channel. The ^y+jets channel has multijet 
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events and T^+jets as main backgrounds. 

In all channels contributions from multijet events are 
modeled using data. Most other background sources 
are modeled through MC simulation. In Run I, both 
collaborations used vecbos [7|| to model VK+jets events. 
VECBOS is a precursor of ALPGEN and provides one of 
the first models of events with many high-momentum 
final-state partons. pythia was used to model Z+jets, 
Drell-Yan, and diboson processes. Background from 
events with single top quark was negligible. In Run II, 
both collaborations used ALPGEN for the simulation of 
the W^+jets background. The treatment of heavy-flavor 
jets is implemented more accurately in ALPGEN and 
parton-jet matching also improves the simulation. For 
the Z+jets background, CDF uses pythia and DO uses 
ALPGEN. For dibosons, both collaborations use PYTHIA. 
Processes with single top-quark are modeled by CDF 
using madevent (76| (based on madgraph \7W). and 
by DO with singletop [zl| (based on comphep [79j]). 

The uncertainty in the description of the VF+jets 
background has three main components: (i) the 
uncertainty on the scale Q 2 , which affects both the overall 
normalization and the differential jet distributions in 
pseudorapidity r\ (80j and pt', (ii) the uncertainty in the 
correction for flavor content of jets to higher order; and 
(iii) the limitation in the MC model we are using to 
reproduce the jet pt and rj distributions in data at low 
Pt and large \rj\. 

E. Jet properties 

After the top quarks decay, the final-state quarks 
and gluons hadronize to produce multiple charged and 
neutral particles that traverse the central tracking 
systems into the calorimeters, where they produce many 
lower-momentum particles through interactions in the 
absorbers of the calorimeters. The observed particles 
tend to cluster in jets that can be assigned to the initial 
partons. For jet reconstruction, the CDF collaboration 
uses a clustering algorithm in (77, 4>) space [8l[ with a cone 
radius of 

CDF K = Vi^v) 2 + (A</>) 2 = 0.4, 

where 4> is the azimuthal angle around the beamline, r\ 
is the pseudorapidity, and A77 or A</> are the widths of 
the cone. DO uses a mid poin t iterative seed-based cone 
algorithm in (y, (f>) space [82J with a radius defined by 

DO Tl = V(Ay) 2 + (A0) 2 = 0.5, 

where the rapidity y = l/2\n((E + p L ) / (E~p L )), E 
is the jet energy, and ph is its longitudinal momentum 
component. 

The jet energy resolution in the central region 
(|ry| < 1) is approximately the same for CDF and 
DO; for CDF it is a(E T )/E T = hO%/yjE T {GeV) © 
3%. For jets in the forward region, however, the 
energy resolution at DO is similar to that in the central 
region, while at CDF it is not as good (ct(Et)/ Et — 



70%/ ^E T (GeV) © 4%). CDF's calorimeter covers 
\r)\ < 3.8 whereas DO's calorimeter covers |?7| < 4.2. 
The DO calorimeter is more homogeneous, so that the 
imbalance in transverse momentum (see Section III G[) 
usually has better resolution at DO. For both CDF and 
DO, to reject jets with mismeasured energy, selections 
on energy deposition are required when clustering the 
energy from the calorimeter cells into jets. When a 
muon is reconstructed within the jet cone, a correction is 
applied to the jet energy to account for the muon and its 
associated neutrino assumed to arise from heavy-quark 
decay. 

Jet energy scale calibrations are applied after 
jet reconstruction. CDF calibrates the transverse 
momentum using test-beam data and single-particle 
simulated events and corrects the jet energy to the parton 
level. Consequently, CDF docs not calibrate the jet 
energy scale in MC events. DO calibrates the energy using 
photon+jets and two-jet data and calibrates jets in data 
as well as in MC to the observed particle level. Particle 
jets arc clustered from stable particles after fragmen- 
tation, including particles from the underlying event, but 
excluding undetected energy from muons and neutrinos. 

CDF's jet calibration [83[ applies two scale factors 
and three offsets to convert the measured transverse 
momentum of a jet to that of the parton that initiated 
the jet. DO's jet calibration [84[ applies three scale 
factors and one offset to the jet energy to convert to the 
particle jet energy scale. The calibrations are expressed 
as follows: 

r<TW parton Pt ^rel — /-< , r~< 

CDF p T = -± Cue + Cbc, 

-ftabs 

D q ^particle _ E3Ct ~ CMLUE 
-Rabs Riel Foe 

The absolute response i? a bs corrects for energy lost in 
uninstrumcntcd regions between calorimeter modules, for 
differences between electromagnetically and hadronically 
interacting particles, as well as for module-to-module 
irregularities. The relative response i? re i is a scale factor 
that corrects forward relative to central jets and Cmi 
is a correction for multiple interactions in the same 
bunch crossing. The function Cue is a correction for 
the jet energy added from the underlying event. DO 
has one offset correction, Cmi,ue, which includes the 
effects of multiple interactions, the underlying event, 
noise from radioactive decays of the uranium absorber, 
and the effect of collisions from previous bunch crossings 
(pile up). The functions Coc and Foe are corrections 
for shower particles scattered in or out of the cone of 
radius 1Z. CDF's correction accounts for MC modeling 
that affects how the parton energy is translated into 
particle jet energy, whereas DO's correction accounts 
for a detector effect caused by the finite cell size in 
the calorimeter coupled with the cone size for the jet 
algorithm. The combined jet energy scale corrections 
increase the measured jet energies by about 20-50%, 
depending on pt and rj. 
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The overall uncertainties on the jet energy scale 
corrections vary from about 2.7% for CDF and 1.1% for 
DO for central jets of transverse energy of 100 GcV to 
3.3% for CDF and 2.2% for DO for forward jets. Central 
jets of 25 GeV have correction uncertainties of 5.9% 
for CDF and 1.4% for DO. For both experiments, the 
uncertainty on the corrections for absolute response i? a bs 
dominate these uncertainties. 

At DO, the jet energy resolution in data is inferior 
than predicted by the detector simulation. Therefore, 
the energies of MC jets are smeared so that the resulting 
resolution in MC matches that in data. Similarly, the 
reconstruction efficiency for jets in data is lower than is 
predicted by the detector simulation, so an appropriate 
fraction of MC jets are randomly removed. Both effects 
are corrected for as functions of jet p? and pseudo- 
rapidity. 

DO Run II analyses include an energy correction to 
simulated jets that depends on jet flavor. There are 
corrections for b jets, other-quark flavor jets (it, d, 
s, and c), and gluon jets implemented in both the 
lepton+jets and dilepton analyses. Such corrections 
refine the simulation by improving the matching of jet 
energies in MC to data. The differences arise from the 
varying electromagnetic fractions and widths of the jets. 
The corrections depend on jet transverse energy and 
pseudorapidity and range from —6% to +2% (27| . 

Both collaborations perform an in-situjet energy scale 
calibration in lepton+jets events for the matrix-element 
mass extraction of mt, and in CDF's alljets and ^T+jets 
measurements of mt ■ The invariant mass of the two jets is 
constrained to a Breit-Wigner distribution for the W 
qq' decay, set to the world average value for the IV-boson 
mass The energies of all jets in the event are then 
rescaled to complete this calibration. 

F. &-quark jet properties 

To separate top-quark events from background and to 
decrease the ambiguity in jet-to-parton matching, it is 
important to identify 6-quark jets. Every ti event has 
two b jets, whereas such jets are rare in background. As 
B hadrons have a mean lifetime of « 10 -12 seconds, 
b jets can be tagged through secondary vertices of 
the B decay a few mm away from the primary pp 
interaction. CDF's 6- tagging algorithm uses the signif- 
icance of the displacement of the secondary vertex in the 
transverse (r 4>) plane for the lepton+jets and ^T+jets 
channels [59|, as well as a jet-probability algorithm for 
|?T+jcts events (63|. One parameter defines the signif- 
icance of the separation of the primary and secondary 
vertices for events with one and two b jets. For jets 
that are within the acceptance of the silicon microstrip 
tracker (i.e., "taggable" jets), this algorithm identifies 
50% of real b jets and 9% of real charm jets, while falsely 
tagging 1% of light jets. DO tags jets by combining 
nine track and secondary-vertex-related variables using 
a neural network [62| . For jets within the acceptance 
of the silicon microstrip detector, this yields efficiencies 



of 65% and 20% for real b and charm jets, respectively, 
while falsely tagging 3% of light jets. 

To identify heavy flavor jets in data and in MC events, 
the tagging algorithm is applied by CDF and DO directly 
to the jets, except for simulated V^+light jets events, 
where CDF uses tag-rate functions measured in multijet 
data, since the rate for directly-tagged MC events is very 
low. After applying direct tagging to b and c jets in 
MC events, DO corrects the tagging efficiencies to match 
those observed in data by randomly dropping the tagging 
of 13% of such jets. For light-flavor jets, DO assigns a a 
per jet mistag weight. 

G. Properties of other event observables 

The uncertainty on m t depends not only on an 
accurate measurement of jet energies and proper 
assignment of flavor, but also on the reconstruction and 
calibration of the other elements of the event, including 
electrons, muons, and the imbalance in transverse 
momentum, taking into account the presence of any 
simultaneous pp interactions in the same bunch crossing. 

The mean number of pp collisions per bunch crossing 
is w 2 in Run I and k 5 in Run II. Such additional 
collisions affect the observed characteristics of the hard 
scatter of interest, and must be included in the MC 
simulation. These extra collisions result mostly in 
the production of low-p^ particles. CDF simulates 
such additional interactions using the pythia model of 
minimum-bias events and overlays them onto the hard 
scatters using a Poisson mean appropriate to the instan- 
taneous luminosity of the data. In a similar manner DO 
overlays randomly-triggered data events with the same 
luminosity profile as the data onto the MC simulated 
events. 

Electrons are identified by matching clusters of 
energy deposited in the electromagnetic layers of the 
calorimeters with tracks that point from the primary 
collision vertex to the clusters. The spatial shapes of the 
showers must agree with those expected for electrons, as 
studied in test beam data. The energy of an electron is 
determined as a combination of the total energy of the 
cluster and the momentum measured from the curvature 
of the matching track. The reconstruction efficiency is 
determined using Z —} ee data by identifying one tight 
charged lepton as tag and using the other charged lepton 
as a probe (tag-and-probe method). The electron energy 
is also recalibrated using such Z events. 

Muons are reconstructed from a central track and 
matched to a track in the outer muon chambers. In 
DO, both the inner and outer trajectories pass through 
magnetic fields and so the transverse momenta of the 
two are therefore required to match. The reconstruction 
efficiency and calibration of px are determined using 
a tag-and-probe method applied on J/ip — > fifj, and 
Z — > pp, events in a manner similar to that used for 
electrons. 

As indicated above, all ti decay channels except for 
alljets events have a large $t- All jet energy calibration 
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corrections arc also propagated to $t in each event. 



III. COMBINATION OF MASS 
MEASUREMENTS 

A. BLUE combination method 

The basic idea of the technique, called the best linear 
unbiased estimator (BLUE) method [34T - [36l ] . used to 
obtain the combined mass m™ mb , an "estimator" of the 
true mass m* ruo , is to calculate a linear weighted sum of 
the results from separate measurements: 



12 
i=l 



(1) 



The m\ are the twelve CDF and DO measurements i of 
mt and 



12 

i=l 



(2) 



The weights are determined using the value of m c t omb that 



true value m* rue 



(m™ mb -m* rue )' 



Variance (m c t omh ) + [Bias(m™ mb )] , (3) 

where the two terms represent the weighted variance and 
bias in the twelve input m t values with 

12 

Variance (mt° mb ) w} Variance(mj) , (4) 



i=l 



and 



Variance(mJ) = [c(tiJ)]' 



(5) 



where a (ml) are the uncertainties on the twelve input 
values given in Table HI 

On average, we expect the input mass measurements 
to be unbiased and we therefore assume 



Bias(m™ mb ) = J2 w i Bias (to j 



(6) 



Equation (3) shows that the BLUE method defines the 
best estimate through a minimization of the variance of 
TOt for an assumed unbiased set of measurements. The 
minimum corresponds to setting the weights to 



W; 



1 /Variance (to J ) 
J2i=i 1 /Variance (toJ) 



(7) 



for uncorrelated input values. Since the input m t values 
arc correlated, the variance in Eq. (4) has to be replaced 
with a covariancc matrix: 

Variance (m 4 comb ) = 

12 12 

u>i Wj Covariance^TOj, , (8) 

i=i j=i 
which is defined as 



Covariance! mj, m\ 



Minimizing Eq. (3) yields 

yV —1 Covariance^ 1 (m\, m\ \ 

W i = 7 A ; 

Ei=i Y,]ti Covariance" 1 [m\, m 3 t j 



(9) 



(10) 



minimizes the squared difference relative to the unknown matrix) , and 



where Covariance \mlml) arc the elements of the 
inverse of the covariance matrix (also known as the error 



Covariance 



a(mfy(ml) (11) 



with Correlation 



the correlation coefficient 



between m\ and m\ . The following sections show 
how the correlation matrix is derived by examining 
the uncertainty components and their individual 
correlations. 

B. Measurement uncertainties 

The uncertainty on any m t measurement has a 
statistical component from the limited number of 
events available for the measurement and a systematic 
component from the uncertainties assigned to the 
calibration of input quantities, to the model of the signal, 
and to the calibration of the mass extraction method. 
Since the first measurements of mt [1, @ , the systematic 
component has been slightly larger than the statistical 
one. As more data became available, the statistical 
uncertainties on m t improved as did the calibrations 
of systematic uncertainty, and the two components 
therefore improved together. 

The systematic uncertainty on each m t measurement 
in this combination is divided into 14 parts. Some 
of them have origin in only one source whereas others 
include several related sources of uncertainties. For 
the latter the patterns of correlation among different 
channels, Tcvatron Run I and Run II, or experiments 
are the same for all sources included in these systematic 
components. The uncertainty on jet energy scale (JES), 
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on the other hand, is split into seven components, which 
do not appply to all measurements, given the significantly 
different approaches to jet energy calibration between 
CDF and DO and the change in the DO procedure between 
Run I and Run II. 

Table [TT] gives the uncertainty of each of the twelve 
top quark mass measurements for the different contri- 
butions to uncertainty, and their effect on the final 
combination. The components of uncertainty are 
defined in the following and can be classified as 
uncertainties in detector response (jet energy scale, 
jet and lepton modeling), uncertainties from modeling 
signal and background (signal modeling, multiple 
interactions model, background estimated from theory 
and background based on data), uncertainties from 
method of mass extraction and statistical uncertainties. 
A detailed description of the methods to evaluate these 
systematic uncertainties is presented in the Appendix. 

1. Jet energy scale 

1.1 Light-jet response (1) 

One subcomponent of the uncertainty in JES 
covers the absolute calibration for CDF's Run I and 
Run II measurements. It also includes small contri- 
butions from the uncertainties associated with modeling 
multiple interactions within a single bunch crossing and 
corrections for the underlying event. 

1.2 Light-jet response (2) 

Another subcomponent of this uncertainty includes 
DO's Run I and Run II calibrations of absolute response 
(energy dependent), the relative response (^-dependent), 
and the out-of-cone showering correction which is 
a detector effect. This uncertainty term for CDF 
includes only the small relative response calibration (77- 
dependent) for Run I and Run II. 

1.3 Out-of-cone correction 

This subcomponent of the JES uncertainty quantifies 
the out-of-cone showering corrections to the MC showers 
for all of CDF's and for DO's Run I measurements that is 
obtained by varying the model for light-quark fragmen- 
tation. 

1.4 Offset 

This subcomponent originates from the offset in DO's 
Run I calibration, which corrects for noise from uranium 
decay, pile-up from previous collisions, and for multiple 
interactions and the model for the underlying event. In 
Run I, the uncertainties are large, but in Run II, owing to 
the smaller integration time for calorimeter electronics, 
they are negligible. CDF's calorimeter does not have the 
same sources of noise and sensitivity to pile-up as DO so 
CDF measurements do not have this term. 

1.5 Model for b jets 

This subcomponent comes from the uncertainty on 
the semilcptonic branching fraction in b decays and from 
differences between two models of 6-jet hadronization. 

1.6 Response to b/q/g jets 



This subcomponent accounts for the difference in 
electromagnetic versus hadronic response of b jets, light- 
quark and gluon jets. CDF corrects for jet flavor as part 
of the main calibration, and defines the uncertainty based 
on the remaining difference in response between b jets and 
light-flavor jets, whereas DO corrects the response for b, 
light-quark (u, d, s, and c) and gluon jets as a function 
of jet pt and r\. 

1.7 In-situ light-jet calibration 

The last part of the uncertainty in jet energy scale is 
from the in-situ calibration of mt- It corresponds to the 
statistical uncertainty from the limited number of events 
used in the fit when using the W-boson mass to constrain 
the energies of the light quarks from the W decay. 

2. Jet modeling 

The uncertainty in jet modeling has two components 
for DO. This uncertainty is negligible for CDF. 

(i) The jet energy resolution is smeared for MC jets 
to match the resolution observed in data and the 
uncertainty on the smearing functions is propagated to 
m t . 

(ii) The identification efficiency in MC events is corrected 
to match that found in data and the uncertainty on the 
correction functions is propagated to m*. 

3. Lepton modeling 

This uncertainty has two components: 

(i) The electron and muon pt scales are calibrated to 
the J/tp and Z-boson mass by both CDF and DO. 
This uncertainty on the calibration is included in the 
measurements of m t . 

(ii) DO smears the muon momentum resolution in MC 
events to match that in data and the uncertainty on this 
correction is included in this term. The uncertainty on 
the electron resolution has a negligible impact on the 
measurements of mt- 

4. Signal modeling 

There are six components to this uncertainty. They 
are combined into one term because the correlations 
between channels arc similar for each component: 

(i) Knowledge of the PDF parametrization. 

(ii) The quark annihilation and gluon fusion fractions 
that differ significantly between leading-log and next-to- 
leading order (NLO) QCD calculations (Run II). 

(iii) The amount of initial and final-state radiation in MC 
signal events differs from that in data and is adjusted 
through the value of Aqcd used in the shower and the 
scales of time and space- like showers. 

(iv) Higher-order QCD corrections to initial and final- 
state radiation differ from precise parton-level models 
and this is not accounted for by the choice of scale for 
the calculations (Run II). 

(v) Our model for jet hadronization is based on angular 
ordering in PYTHIA with Tune A underlying event tuning. 
Parton showering and the underlying event can also be 
simulated with HERWIG and JIMMY (8f|[86[. The effect of 
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TABLE II: The uncertainty in GeV from each component for the twelve measurements of mt and the resulting Tevatron 
combination. The total uncertainties are obtained by adding the components in quadrature. The entries "n/a" stand for 
"not applicable" and "n/e" for "not evaluated." The non-evaluated uncertainties were not considered as significant sources of 
uncertainty for Run I measurements. 




Channel 


Run 


Exp. 




Jet 


energy scale systematics 








Other systematics 














Lepton+jets 


II 


CDF 


0.41 


0.01 


0.27 


n/a 


0.23 


0.13 


0.58 


0.00 


0.14 


0.56 


0.10 


0.27 


0.06 


0.10 


0.65 


0.80 


0.67 


1.23 


Lcpton+jets 


11 


DO 


n/a 


0.63 


n/a 


n/a 


0.07 


0.26 


0.46 


0.36 


0.18 


0.77 


0.05 


0.19 


0.23 


0.16 


0.83 


0.83 


0.94 


1.50 


Lcpton+jets 


I 


CDF 


3.4 


0.7 


2.7 


n / a 


0.6 


n/e 


n/a 


n/c 


n / c 


2.7 


n / c 


1.:! 


n/e 


0.0 


5.1 


4.4 


2.8 


7.3 


Lcpton+jets 


I 


DO 


n/a 


2.5 


2.0 


1.3 


0.7 


n/e 


n/a 


n/c 


n/e 


1.3 


n/e 


1.0 


n/e 


0.6 


3.6 


3.5 


1.6 


5.3 


Alljets 


II 


CDF 


0.38 


0.04 


0.24 


n/a 


0.15 


0.03 


0.95 


0.00 


n/a 


0.64 


0.08 


0.00 


0.56 


0.38 


1.43 


1.06 


0.91 


2.00 


Alljcts 


I 


CDF 


4.0 


0.3 


3.0 


n/a 


0.6 


n/e 


n/a 


n/c 


n/a 


2.1 


n / c 


1.7 


n/e 


0.6 


10.0 


5.0 


2.6 


11.5 


Dilcptons 


11 


CDF 


2.01 


0.58 


2.13 


n/a 


0.33 


0.14 


n/a 


0.00 


0.27 


0.80 


0.23 


0.24 


0.14 


0.12 


1.95 


3.01 


0.88 


3.69 


Dilcptons 


II 


DO 


n/a 


0.56 


n/a 


n/a 


0.20 


0.40 


0.55 


0.50 


0.35 


0.86 


0.00 


0.00 


0.20 


0.51 


2.36 


0.90 


1.11 


2.76 


Dilcptons 


I 


CDF 


2.7 


0.6 


2.6 


n/a 


0.8 


n/e 


n/a 


n/c 


n/c 


3.0 


n / c 


0.3 


n/e 


0.7 


10.3 


3.9 


3.0 


11.4 


Dilcptons 


I 


DO 


n/a 


1.1 


2.0 


1.3 


0.7 


n / e 


n/a 


n/c 


n / c 


1.9 


n/e 


1.1 


n/e 


1.1 


12.3 


2.7 


2.3 


12.8 


#T+jets 


II 


CDF 


0.45 


0.05 


0.20 


n/a 


0.00 


0.12 


1.54 


0.00 


n/a 


0.78 


0.16 


0.00 


0.12 


0.14 


1.80 


1.64 


0.78 


2.56 


Decay length 


11 


CDF 


0.24 


0.06 


n/a 


n/a 


0.15 


n/e 


n/a 


0.00 


n / a 


0.90 


0.00 


0.80 


0.20 


2.50 


9.00 


0.25 


2.80 


9.43 


Tevatron Combination 


0.12 


0.19 


0.04 


0.00 


0.15 


0.12 


0.39 


0.11 


0.10 


0.51 


0.00 


0.14 


0.11 


0.09 


0.56 


0.49 


0.57 


0.94 



the difference on mt between the two models is included 
in this term. 

(vi) Final-state partons and remnants of the protons and 
antiprotons are connected through color strings, which 
affect the distributions of jets. Since this effect is not 
included in the model for tt signal, the value of mt has 
an uncertainty from this omission (Run II). 

5. Multiple interactions model 

The number of soft pp events overlaid on each MC 
event has a Poisson distribution. The mean number 
does not equal exactly the number seen in data since the 
luminosity increased as the Tevatron run progressed. The 
top-quark mass is measured as a function of the number 
of multiple interactions in signal events by CDF and the 
signal MC events are reweighted to match the distri- 
bution seen in data by DO and the related uncertainties 
are included here. 

6. Background from theory 

There are four components in this uncertainty: 

(i) Difference between NLO calculations of the fraction of 
heavy-flavor jets in W^+jets events. The ALPGEN model 
underestimates this fraction. 

(ii) Impact of factorization and rcnormalization scales 
on the VF+jets simulation, which affects the background 
model for distributions characterizing jets. 

(iii) The theoretical cross sections used to normalize all 
MC estimated background processes (except for IF+jcts 



for CDF and DO lepton+jets measurements, and Drell- 
Yan production for CDF dilepton measurements), 
(iv) Impact of difference between the MC modeling of 
background kinematic distributions and those observed 
in data. 

7. Background based on data 

This refers primarily to uncertainties from the normal- 
ization of certain background components to data. These 
include multijet backgrounds in the lepton+jets, alljcts, 
and ^T+jets analyses, the W^+jets background in the DO 
lcpton+jets analyses, and the Drell-Yan backgrounds in 
the CDF dilepton analyses. 

DO also considers the following four components of 
uncertainty: 

(i) the uncertainty from correcting the MC events to 
match the trigger efficiency in data which is based on 
the turn-on response for each trigger element. 

(ii) the uncertainty from applying tag-rate and 
taggability corrections to MC events to make the 
efficiencies match the data for each jet flavor. 

(iii) the uncertainty on the tt signal fractions in the 
samples used to calibrate the measurements. 

(iv) the uncertainty on the fraction of multijet events 
included in the pseudocxperiments used for calibration. 

8. Calibration method 

The extracted values of m t are calibrated using a 
straight-line fit to the relationship between input mass 



TABLE III: Correlations in systematic uncertainties (in percent) among the different measurements 
of nit. 
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TABLE IV: Correlations in systematic uncertainties (in percent) among the different measurements 
of mt (continued). 
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and measured mass in simulated pseudoexperiments. 
This term includes the systematic uncertainties from the 
slope and offset of this calibration. 

9. Statistical uncertainty 

The statistical uncertainties are determined from 
the number of data events in each of the twelve 
measurements. 

Figure [5] shows the relative contribution for each major 
uncertainty to the analysis channels in Run II. The 
Appendix provides more detail on how each of the sources 
of the uncertainties is estimated. 
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FIG. 5: The average uncertainties for CDF and DO for 
each Run II measurement and for the Tevatron combination, 
separated according to major components (see Table ["VTTTl in 
the Appendix for details on the systematic categories. In this 
figure, the jet and lepton modeling systematic uncertainties 
are grouped into the modeling background category.). 



C. Uncertainty correlations 

Tables IIIII and IIVI indicate how uncertainties are 
correlated between measurements. There are seven 
patterns of correlation: 

(i) Statistical uncertainty and calibration method 
uncertainty are not correlated among the measurements. 

(ii) Correlations among DO measurements that 
implement the same final jet energy corrections for 
the uncertainty from in-situ light-jet calibration. 

(iii) Correlations among CDF measurements that use the 
same data samples for the uncertainty from background 
based on data. 

(iv) Correlations among all measurements in the same 
tt decay channel for the uncertainty from background 
estimated from theory. 

(v) Correlations of measurements within the same 
experiment for a given run period for the uncertainties 
from light-jet response (2), offset, response to b/q/g jets, 
jet modeling, lepton modeling and multiple interactions 
model. 



(vi) Correlations for measurements within the same 
experiment such as the uncertainty from light-jet 
response (1). 

(vii) Correlations among all measurements such as the 
uncertainties from out-of-cone correction, model for b jets 
and signal modeling. 

We assume that all sources correspond to either no 
or 100% correlation. A check of this assumption (see 
Section I IV Bp shows that it has negligible effect on the 
combined value and uncertainty of mt- 

D. Measurement correlations 

The uncertainties shown in Table [TT] and their 
correlations shown in Tables IIIII and IIVI provide the 
correlations among the twelve input values of m t . The 
correlation matrix for these measurements, as returned 
by the combination procedure, is shown in Table [V] 
The inversion of the covariance matrix built with the 
correlation matrix defines the measurement weights, as 
described in Section lillAI 

E. Measurement weights 

As discussed in Section MI Al the combined mass 
TO comb j g defined through the set of weights that minimize 
the squared difference between r?"i£ omb and the true value 
of m t , which is equivalent to minimizing the sum of the 
covariance matrix elements. Table [V] gives the weights 
u>i for each of the input measurements as determined 
in this minimization. A weight of zero means that an 
input measurement has no effect on m™ b . The Run I 
measurement weights are negative which reflects the fact 
that the correlations for these and other measurements 
are larger than the ratio of their total uncertainties [351 ]. 
In this case, the less precise measurement may acquire 
a negative weight. Input measurements with negative 
weights still affect the value of m£ omb and reduce the 
total uncertainty. By design, the sum of the weights is 
set to unity. 



IV. RESULTS OF THE COMBINATION 

A. Tevatron top-quark mass result 

Combining the twelve independent measurements of 
mt from the CDF and DO collaborations yields 



n comb 



173.18 ± 0.56 (stat) ±0.75 (syst) GeV 
173.18 ±0.94 GeV. 



The uncertainties are split into their components in 
Table [TTJ and Fig. [5j The jet energy scale contributes 
0.49 GeV to the total systematic uncertainty. Of this, 
0.39 GeV arises from limited statistics of the in-situ 
JES calibration and 0.30 GeV from the remaining contri- 
butions. Figure |H] summarizes the input m t values and 
the combined result. 
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TABLE V: Correlations in % among the input mt measurements and their weights in the BLUE combination. 
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FIG. 6: The twelve input measurements of mt from the Tevatron collider experiments along with the resulting combined value 
of m£ omb . The grey region corresponds to ±0.94 GeV. 



We assess the consistency of the input mt test statistic, defined as follows: 
measurements with their combination using a \ 2 

Xcomb = K - ™r mb f 



x Covariance 1 {m\ 1 ml S j (m{ - m£ omb ^ , 
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TABLE VI: Separate calculations of raf 
experiment, and their x 2 probabilities. 



for each tt decay mode, by run period, and by 



Subset 



Consistency \ 
(Degrees of freedom — 1) 



X" probability 



h 
"-4 



En 
Q 



C 



Lepton+jets 173.4 ± 1.0 



Alljets 

Dileptons 

^T+jets 

Run II 

Run I 

CDF 

DO 



172.7 ± 1.9 
171.1 ± 2.1 
172.1 ± 2.5 
173.6 ± 1.0 
180.0 ±4.1 
172.5 ± 1.0 
174.9 ± 1.4 



0.14 
1.51 
0.28 



0.40 
0.04 



1.51 
0.40 



0.12 



0.28 
0.04 
0.12 



2.56 



71% 22% 60% 

71% 53% 85% 

22% 53% 73% 

60% 85% 73% 



<)% 



11% 



where ml is a column vector of the twelve m t 
inputs, m£ omb is a matching column vector for the 
measurements adjusted in the previous minimization, 
and the superscript T denotes the transpose. We find 

Xcomb = 8-3 for 11 degrees of freedom, 

which is equivalent to a 69% probability for agreement 
(i.e., p- value for the observed x 2 value) among the twelve 
input measurements. 



B. Consistency checks 

We check one aspect of the assumption that biases 
in the input m t are on average zero (see section IIII Ap 
by calculating separately the combined m™ mb for each 
tt decay mode, each run period, and each experiment. 
The results are shown in Table I VII The resulting m™ mb 
values are calculated using all twelve input measurements 
and their correlations. The \ 2 test statistic provides 
the compatibility of each subset with the others and is 
defined as: 

2 

Xsubl,sub2 

(mj ubl - m^ ub2 ) 2 Covariance" 1 (m™ bl , mf lb2 ) . 

The x 2 values in Table I VII show that biases in the input 
measurements are not large. 

To check the impact of the assumption that the 
systematic uncertainty terms are either 0% or 100% 
correlated between input measurements, we change all 
off-diagonal 100% values to 50% (see Tables IIIII and 
IIVP and recalculate the combined top-quark mass. This 
extreme change shifts the central mass value up by 
0.17 GeV and reduces the uncertainty negligibly. The 
chosen approach is therefore conservative. 



C. Summary 

We have combined twelve measurements of the mass 
of the top quark by the CDF and DO collaborations at 
the Tevatron collider and find: 

m comb = 173 18 ± 5g ( gtat ) ± Q 75 ( gyst ) Ge y 

which corresponds to a precision of 0.54%. The result 
is shown in Table IVHI together with previous combined 
results for comparison. The input measurements for 
this combination use up to 5.8 fb _1 of integrated 
luminosity for each experiment, while 10 fb _1 are now 
available. We therefore expect the final combination 
to improve in precision with the use of all the data, 
but also from analyzing all tt decay channels in both 
experiments, and from the application of improved 
measurement techniques, signal and background models, 
and calibration corrections to all channels that will 
reduce systematic uncertainties. Currently, there are also 
some overlaps of the systematic effects that are included 
in different uncertainty categories. In addition to the in- 
still light-jet calibration systematic uncertainty that will 
scale down with the increase of analyzed luminosity, these 
levels of double counting are expected to be reduced for 
the next combination. The combination presented here 
has a 0.54% precision on m* making the top quark the 
particle with the best known mass in the SM. 
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TABLE VII: Mass measurements of the top quark from 1999 
until this publication at the Tevatron collider. 
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APPENDIX: EVALUATION OF SYSTEMATIC 
UNCERTAINTIES 

Systematic uncertainties arise from inadequate 
modeling of signal and backgrounds, and the inability to 
reproduce the detector response with simulated events. 
Systematic uncertainties also arise from ambiguities in 
reconstructing the top quarks from their jet and lepton 
remnants. We minimize such uncertainties by using 
independent data to calibrate the absolute response of 
the detector, and we use state-of-the-art input from 
theory for modeling the signal and backgrounds. We use 
alternative models for signal and different parameters 
for modeling backgrounds to check our assumptions. 

Tabic IVIIII lists the uncertainties from the Run II 
lepton+jets measurements for CDF and DO that arc 
based on the matrix element technique (2(| |27|- These 
two measurements provide most of the sensitivity to the 
combined mt result and are discussed below. Before 
explaining how each individual systematic uncertainty is 
estimated, we will first discuss how the uncertainties from 
different sources are propagated to mt and how they are 
calculated using ensembles of pscudoexperimcnts. 

Uncertainties related to the performance of the 
detector and calibration of the reconstructed objects, 
such as jet energy scale (JES), the modeling of jets, 
leptons, and triggers, and calibration of the ^-tagging 
algorithms, are evaluated by shifting the central values 
of their respective parameters by ±1 standard deviations 



(a) that correspond to the uncertainties on each value. 
This is done using Monte Carlo (MC) ti events for nit = 
172.5 GeV. The integrations over the matrix element 
are performed again for each shifted sample and define 
shifts in mt that correspond to each independent source 
of systematic uncertainty. These uncertainties are not 
determined at other m t values and it is assumed that 
their dependence on mt is minimal. 

For uncertainties that arise from ambiguities in the 
modeling of the ti signal, which include the uncertainties 
from initial and final-state radiation, higher-order QCD 
corrections, 6-jet hadronization, light-jet hadronization, 
the underlying event model, and color reconncction, we 
generate simulated ti events using alternative models 
also at m t = 172.5 GeV. These events are processed 
through detector simulation, arc reconstructed and the 
probability density is calculated by integration over the 
matrix elements. 

For the uncertainties from the choice of parton distri- 
bution functions, the ratio of contribution from quark 
annihilation and gluon fusion, and models for overlapping 
interactions, we rewcight the fully reconstructed 
simulated ti MC events at m t = 165, 170, 172.5, 175, and 
180 GeV to reflect the uncertainty on ±1ct range on each 
parameter, and extract its impact on mt- 

Each method used to measure m t is calibrated 
using ti MC events generated at m t = 
165,170,172.5,175,180 GeV, which provide the 
relationship between input and "measured" masses. 
A straight line is fitted to these values, representing 
a response function that is used to correct the m t 
measurement in data. 

Systematic uncertainties are evaluated using studies 
of ensembles of pscudoexperiments. For each of the 
shifted or rcweightcd sets of events, and those based on 
alternative models or different generated m t , we create an 
ensemble of at least 1000 pscudoexperiments, by means 
of binomially-smeared signal and background fractions 
that match the expectation in the data sample and with 
the total number of events in each pseudoexperiment 
equal to the number of events observed in data. We use 
the ensembles of such pseudoexperiments to assess the 
difference between generated and measured mass, and to 
calibrate the method of mass extraction. 

For the uncertainty on background, we change the 
fraction of background events in the pseudoexperiments 
within their uncertainties and remeasure the top-quark 
mass. 

For the BLUE combination method, the uncertainties 
must be defined symmetrically around the central mass 
value, and this requirement determines part of the 
following definitions of uncertainty. 

For the uncertainties obtained in ensemble studies with 
shifted or reweighted parameters, m^ corresponds to the 
+1<7 shift in the input parameter and m^ corresponds to 
the —1(7 shift. The systematic uncertainty on the value of 
m t from these parameters is defined as ± \mf — to^I/2, 
unless both shifts arc in the same direction relative to the 
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TABLE VIII: Individual components of uncertainty on CDF and DO mt 
measurements in the lepton+jets channel for Run II data [26l . |27| . 





Uncertainty [G 1 


sV] 


Systematic 


CDF (5.6 fb" 1 ) 


DO 


(3.6 fb- 1 ) 


Source 


m t = 173.00 GeV 


m t = 


174.94 GeV 



DETECTOR RESPONSE 
Jet energy scale 



Light-jet response (1) 


n A i 


n/a 


Light- j ct response (2) 


n ni 

U.Ul 


U.DO 


Out- of- cone correction 


27 




iVloclcl for b jets 


0.23 


0.07 


Semileptonic b decay 


0.16 


0.04 


b-jet hadronization 


0.16 


0.06 


Response to b/q/g jets 


0.13 


0.26 


In-situ light-jet calibration 


0.58 


0.46 


Jet modeling 


0.00 


0.36 


Jet energy resolution 


0.00 


0.24 


Jet identification 


0.00 


0.26 


Lcpton modeling 


0.14 


0.18 


MODELING SIGNAL 






Signal modeling 


0.56 


0.77 


Parton distribution functions 


0.14 


0.24 


Quark annihilation fraction 


0.03 


n/a 


Initial and final-state radiation 


0.15 


0.26 


Higher-order QCD corrections 


n/a 


0.25 


Jet hadronization and underlying event 


0.25 


0.58 


Color reconnection 


37 


28 


Multiple interactions model 


0.10 


0.05 


MODELING BACKGROUND 






Background from theory 


0.27 


0.19 


Higher-order correction for heavy flavor 


0.03 


0.07 


Factorization scale for W-hjets 


0.07 


0.16 


Normalization to predicted cross sections 


0.25 


0.07 


Distribution for background 


0.07 


0.03 


Background based on data 


0.06 


0.23 


Normalization to data 


0.00 


0.06 


Trigger modeling 


0.00 


0.06 


b-tagging modeling 


0.00 


0.10 


Signal fraction for calibration 


n/a 


0.10 


Impact of multiset background on the calibration 


n/a 


0.14 


METHOD OF MASS EXTRACTION 






Calibration method 


0.10 


0.16 


STATISTICAL UNCERTAINTY 
UNCERTAINTY ON JET ENERGY SCALE 
OTHER SYSTEMATIC UNCERTAINTIES 


0.65 
0.80 
0.67 


0.83 
0.83 
0.94 


TOTAL UNCERTAINTY 


1.23 


1.50 



CDF's calibration of the absolute jet energy scale uses 
the single-pion response to calibrate jets in data and 
to tune the model of the calorimeter in the simulation. 
Uncertainties of these processes form the greatest part 
of the JES uncertainty. Small constant terms are 
added to account for the model of jet fragmentation, for 
calorimeter simulation of electromagnetically decaying 
particles, and to take into account small variations of 
the absolute calorimeter response over time. The total 
resulting uncertainty on the absolute JES is 1.8% for 
20 GeV jets rising to 2.5% for 150 GeV jets. 

At high Tevatron instantaneous luminosities, more 
than one pp interaction occurs during the same bunch 
crossing, and the average number of interactions depends 
linearly on instantaneous luminosity and changed from 
w 1 to 8 between the start and the end of Run II. If 
the final-state particles from these extra pp interactions 



nominal value, in which case the systematic uncertainty 
is defined as the larger of \mf — m t \ or \wl[ — m t \. 

For the values obtained from a comparison between 
two or more models, the systematic uncertainty is taken 
as ± of the largest difference among the resulting masses 
(without dividing by two). 

1. Jet energy scale 

The following seven terms (1.1 - 1.7) refer to the jet 
energy scale: 

1.1 Light-jet response (1) 

This uncertainty includes the absolute calibration of 
the CDF JES for Run I and Run II and the smaller effects 
on JES from overlapping interactions and the model for 
the underlying event. 
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overlap with the jets from a tt event, the energy of 
these jets is increased, thereby requiring the correction. 
The uncertainty on this correction depends on vertex- 
reconstruction efficiency and the rate for misidentifying 
vertices. The impact of these effects is checked on data 
samples, including W — > ev 1 minimum bias, and multijet 
events with a trigger threshold of 100 GeV. CDF finds 
an uncertainty of 0.05 GeV per jet. This uncertainty was 
estimated early in Run II. With increasing instantaneous 
luminosity, this correction was insufficient, and another 
systematic uncertainty term was introduced through the 
"multiple-interactions-model" term, which is described 
later. 

CDF includes the impact of the underlying event on 
JES in this component of uncertainty. The proton and 
antiproton remnants of the collision deposit energy in 
the calorimeter, and these can contribute to the energy 
of the jets from tt decay, which must be subtracted before 
m t can be measured accurately. CDF compares the 
"Tune A" underlying event model [73| in pythia 
with the jimmy model [H, [H| in herwig [13 using 
isolated tracks with pt > 0.5 GeV. The data agree well 
with Tune A, which is expected since it was tuned to 
CDF data, but differ from jimmy by about 30%. This 
difference is propagated to the absolute calibration of 
JES and yields a 2% uncertainty for low-pr jets, and 
less than 0.5% for 35 GeV jets. 

MC tt events are generated by CDF with jet energies 
shifted by the above three uncertainties, and the resulting 
shifts in m t are used to estimate the uncertainty. The 
overall uncertainty on m* from these combined sources is 
0.24% for lepton+jets, 0.22% for alljets, 1.18% for CDF 
Run II dilepton data, and 0.26% for ^T+jets for Run II 
data of CDF. 

1.2 Light-jet response (2) 

This uncertainty term represents almost all parts of 
DO Run I and Run II calibrations of JES. The absolute 
energy scale for jets in data is calibrated using 7+jct 
data with photon px > 7 GeV and |f? 7 | < 1.0, and jet 
Pt > 15 GeV and \ru e t\ < 0.4, using the projection 
fraction" method j84j. Simulated samples of 7+jcts and 
Z+jets events are compared to data, and used to correct 
the energy scale for jets in MC events. The JES is 
also corrected as a function of rj for forward jets relative 
to the central jets using 7+jets and dijets data. Out- 
of-cone particle scattering corrections are determined 
with 7+jets data and simulated events, without using 
overlays of underlying events, to avoid double-counting 
of this effect. Templates of deposited energy are formed 
for particles belonging to and not belonging to a jet 
using 23 annular rings around the jet axis for TZ(y, <j>) — 
sj (Ay) 2 + (A0) 2 < 3.5. All of these calibration steps are 
combined and the total uncertainty on JES is calculated 
for light jets and heavy-flavor jets (independent of the 
type of jet). The resulting DO uncertainty on m t for 
Run II lepton+jets events is 0.36% and 0.86% for dilepton 
data. 



This uncertainty term also includes the relative jet 
energy correction as a function of jet 77 for CDF. This 
is measured using dijet data, along with pythia and 
herwig simulations of tt events generated with shifted 
jet energies, and lead to the following uncertainties on 
Run II measurements of m t : 0.01% for lepton+jets, 
0.02% for alljets, 0.34% for dileptons, and 0.03% for 
^T+jets. 

1.3 Out-of-cone corrections 

For all CDF measurements and for DO Run I, this 
uncertainty component accounts for energy lost outside 
the jet reconstruction cone, and uses the difference 
between two models of light-quark and gluon fragmen- 
tation and simulation of the underlying event. DO 
changed the way it measures the out-of-cone uncertainty 
between Run I and Run II, and this uncertainty for DO 
Run II measurements is therefore included in the light-jet 
response (2) term, described previously. 

Energy is lost from the cone of jet reconstruction when 
a quark or gluon is radiated at large angle relative to 
the original parton direction, or when the fragmentation 
shower is wider than the cone, or when low momentum 
particles are bent out of the cone by the axial magnetic 
field of the detector. Energy is gained in the cone from 
initial-state radiation and from remnants of spectator 
partons, called collectively at CDF the underlying event. 
The two models compared by CDF in Run II are PYTHIA 
with Tune A for the underlying event, and herwig with 
the jimmy modeling of the underlying event. For the 
narrow cone size of TZ = 0.4 used in measurements of 
m t , more energy is lost from the cone than gained. The 
correction is measured using pythia dijet events and 
data in the region 0.4 < TZ < 1.3. A small constant 
is added to compensate for energy outside the TZ > 1.3 
region ("splash out"). The correction is largest for jets 
at low transverse momentum: +18% for px = 20 GeV 
jets and < 4% for jets with pt > 70 GeV. A detailed 
description of the method can be found in Ref . [83| . 

The uncertainty on these corrections is measured by 
comparing 7 + jets data to the two simulations. The 
largest difference between either of the models and data 
is taken as the uncertainty (the difference between the 
two models is very small). For jets with pt = 20 GeV, 
the uncertainty on the jet energy scale is 6%, and for 
jets above 70 GeV, it is 1.5%. These translate into 
uncertainties on CDF Run II m t measurements of 0.16% 
for the lepton+jets measurement, 0.14% for alljets, 1.25% 
for dileptons, and 0.12% for ^T+jets. 

1.4 Energy offset 

This uncertainty term is specific to DO Run I 
measurements. It includes the uncertainty arising from 
uranium decays noise in the calorimeter and from the 
correction for multiple interaction to JES. These lead to 
uncertainties in m t of 0.72% for lepton+jets and 0.77% 
for dilepton events. In Run II, the integration time for 
the calorimeter electronics is short, after the upgrade to 
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shorter bunch-crossing time (3.5 fis to 396 ns). This effect 
results in a negligible uncertainty on the offset for DO 
Run II measurements of mt ■ 

1.5 Model for b jets 

(i) Semileptonic b decay 

The uncertainty on the semileptonic branching fraction 
(10.69 ± 0.22) x 10~ 2 (PDG 2007 values) of B hadrons 
affects the value of m t . Both collaborations reweight 
ti events by ± the uncertainty on the central value 
(±2.1%), and take half the resulting mass difference as 
the uncertainty on m t : 0.09% for CDF and 0.03% for DO. 

(ii) b-jet hadronization 

For its nominal m t measurements, CDF uses the 
default pythia model of b-jct fragmentation based on the 
Bowler model HH (r q = 1.0, a = 0.3, b = 0.58), where 
r q is the Bowler fragmentation-function parameter and a 
and b are Lund fragmentation function parameters. DO 
uses a model with these parameters tuned to data from 
ALEPH, DELPHI, and OPAL ^ (r q = 0.897 ± 0.013, 
a = 1.03 ± 0.08, b = 1.31 ± 0.08). To measure the 
uncertainty on these models, CDF compares its m t values 
to those measured with the LEP parameters used by 
DO, and to those from the SLD experiment at SLC [89| 
(r q = 0.980 ± 0.010, a = 1.30 ± 0.09, b = 1.58 ± 0.09). 
DO compares the measured m t with the LEP parameters 
to the one from SLC. The resulting uncertainties on the 
mt extracted from the lepton+jets channel are 0.09% for 
CDF and 0.03% for DO. 

For some analyses, the determination of the 
uncertainties in (i) and (ii) may be affected by statistical 
fluctuations of the MC samples. 

1.6 Response to b/q/g jets 

The calibrations of JES described in the first two 
paragraphs of the Appendix are derived on samples 
dominated by "light quark" and gluon jets and applied to 
all jets. However the calorimeter response to heavy-flavor 
jets differs in that these particles often decay semilepton- 
ically and the b jet will have some energy lost through the 
escaping neutrino. Bottom quark jets can also contain 
an electron that showers in a pattern different than for 
hadronic particles, or the jet may contain a muon that 
does neither produce a shower nor gets absorbed in the 
calorimeter. Bottom jets also differ from light jets in the 
distribution of their shower and particle content. Since 
every ti event contains two b jets, it is important to 
understand their energy calibration after the application 
of the previous overall corrections. 

CDF measures an uncertainty from the difference 
between the 6-jets response and light-flavor jets response 
in Run II. CDF takes sets of MC ti events and 
cluster particles into jets classifying each such particle 
jet as a b jet or a light jet [81(. Single-particle 
response for data and for MC events are applied to 
the formed particle jets to predict the energy measured 



in the calorimeter. A double ratio is calculated: 

(pdata /;p MC )h .^/^data/pMC)^ which ig found tQ 

be 1.010. The uncertainty on m t is measured by 
generating new tt samples with the fo-jet scale shifted by 
this 1% difference, which results in 0.1% uncertainty in 
mt for the lepton+jets measurement. 

For Run II measurements, DO corrects the transverse 
momentum distributions of jets differently in four regions 
of detector pseudorapidity to make the MC response 
match that in data (after the main JES calibration) as a 
function of jet flavor: b jets, light-quark jets (u, d, s, c) 
and gluon jets [13] • The correction functions are shifted 
up and down by their uncertainties and the extracted 
shifts in mt are used to define the resulting uncertainty 
on m t of 0.15% for the lepton+jets measurement and 
0.23% for the dilepton measurement. 

1.7 In-situ light-jet calibration 

In ti events where one or both W bosons decay to qq', 
the world average value of My/ 0] is used to constrain the 
jet energy scale for light-quark jets in-situ [9(J[9l|]. CDF 
and DO perform simultaneous measurements of mt and 
Myy, and fit a linear function to the JES for light-quark 
jets that is applied to all the jets to improve precision of 
m t . 

CDF measures the in-situ rescaling factor indepen- 
dently in their lepton+jets, alljets, and ^T+jets analyses 
and so these terms are uncorrelated. DO applies the 
rescaling derived from their lepton+jets measurement 
to dilepton events, and these uncertainties are therefore 
correlated. 

The uncertainty from the in-situ calibration is 
determined through a two-dimensional minimization of 
a likelihood that is a function of top-quark mass and 
JES. The extracted JES is then shifted relatively to 
its measured central value and a one-dimensional fit is 
performed to the top-quark mass. The difference in 
quadrature between the uncertainty on mt from the first 
and second fits is taken as the uncertainty on mt from the 
in-situ calibration, giving 0.34% for CDF's lepton+jets 
measurement, 0.27% for DO's lepton+jets result, 0.55% 
for CDF's alljets, 0.89% for their ^T+jets measurement, 
and 0.32% for DO's dilepton measurement. 

2. Jet modeling 

Applying jet algorithms to MC events, CDF finds that 
the resulting efficiencies and resolutions closely match 
those in data. The small differences propagated to m t 
lead to a negligible uncertainty of 0.005 GcV, which is 
then ignored. DO proceeds as follows. 

(i) Jet energy resolution 

The modeling of the jet energy resolution is corrected 
in DO to match that in data. The value of mt is then 
remcasured using MC samples with jet energy resolution 
corrections shifted up and down by their uncertainties, 
resulting in an uncertainty on mt of 0.18%. 
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(ii) Jet identification 

DO applies correction functions to MC events to match 
the jet identification efficiency in data. The uncertainty 
on m t is estimated by reducing the corrections by la and 
rcmcasuring the mass in the adjusted MC samples. The 
efficiency can only be shifted down and not up because 
jets can be removed from the simulated events but not 
added. The uncertainty on m t is therefore set to ± the 
single-sided shift and is 0.15%. 

3. Lepton modeling 

(i) Momentum scale for leptons 

In Run II, the electron and muon channels for CDF 
and the muon channels for DO are used to calibrate 
the lepton momentum scales by comparing the invariant 
dilepton mass m tU 2 = \J {En + E n ) 2 - (wi + P12) 2 for 
J/ip — s> It and Z ->• It decays in MC events with 
data. The positions of the resonances observed in the 
mu distributions reflect the absolute momentum scales 
for the leptons. CDF and DO perform a linear fit as a 
function of the mean value of transverse momentum to 
the two mass points (3.0969 GeV and 91.1876 GeV Q), 
assuming that any mismatch is due to an uncertainty 
in the calibration of the magnetic field. DO also fits 
a quadratic relation assuming that the difference in 
scale arises from misalignment of the detector. The 
value of nit is measured using MC ti ensembles without 
rescaling lepton pt and with lepton pt values rescaled 
using these fitted relations. Half of the largest difference 
in extracting mt is taken as its systematic uncertainty 
resulting from lepton pr scale. For muon measurements 
from DO, the largest shift is observed for the linear 
paramctrization. In Run I, this source of uncertainty was 
neglected as it was negligible relative to other sources of 
uncertainty. 

In DO Run II measurement of the M^-boson mass in the 
electron decay channel, it was found that 0.26 radiation 
lengths of material was left out in the GEANT modeling 
of the solenoid [92j . The Z-boson mass peak was used to 
calculate a quadratic correction to the electron energy by 
comparing MC events generated with additional solenoid 
material to data. This correction was then propagated 
to the mt measurement. 

The uncertainties on the m; measurements from the 
lepton momentum scale are 0.08% for CDF lepton+jets 
measurements and 0.10% for DO, and 0.16% for CDF 
dilepton measurements and 0.28% for DO dilepton 
results. 

(ii) Lepton momentum resolution 

The muon momenta in simulated events at DO are 
smeared to match the resolution in data. The uncertainty 
on this correction corresponds to an uncertainty on mt 
of 0.17%. 

4. Signal modeling 



(i) Parton distribution functions 

In Run I, the uncertainties from choice of parton distri- 
bution functions (PDF) are determined by measuring 
the chan ge i n m t using the MRS A' set [93f instead of 
MRSDg [6l or CTEQ4M 



and are found to be 
negligible. 

In Run II, the uncertainty is measured by CDF 
by comparing CTEQ5L results with MRST98L 
by changing the value of a s in the MRST98L model, 
and by varying the 20 eigenvectors in CTEQ6M 
The total uncertainty is obtained by combining these 
sources in quadrature. DO measures this uncertainty 
by reweighting the pythia model to match possible 
excursions in the parameters represented by the 20 
CTEQ6M uncertainties, and taking the quadratic sum 
of the differences. The resulting uncertainty on m t is 
0.08% for CDF and 0.14% for DO. 

(ii) Fractional contributions from quark annihilation and 
gluon fusion 

In Run I, this source of uncertainty in ti production 
is not considered. In Run II, CDF estimates the effect 
on m t by reweighting the gluon fusion fraction in the 
PYTHIA model from 5% to 20% [H[ . The uncertainty on 
mt is found to be 0.02%. This uncertainty is included by 
DO in the systematic component (iv) below, where the 
effects of higher-order QCD corrections are discussed. 

(iii) Initial and final-state radiation 

Initial and final-state radiation refers to additional 
gluons radiated from the incoming or outgoing partons 
or from the top quarks. Jets initiated by these gluons 
affect the measured value of m t because they can be 
misidentified as jets from the final-state partons in top- 
quark decay. Extensive checks were performed in Run I 
measurements to assess the effects of initial and final- 
state radiation by varying parameters in herwig. 

In Run II, uncertainties from initial and final-state 
radiation are assessed by both collaborations using a 
CDF measurement [96[ in Drell-Yan dilepton events that 
have the same qq initial state as most ti events, but 
no final-state radiation. The mean pt of the produced 
dilepton pairs is measured as a function of the dilepton 
invariant mass, and the values of Aqcd and the Q 2 scale 
in the MC that matches best the data when extrapolated 
to the ti mass region are found. CDF's best-fit values are 
Aqcd (5 flavors) = 292 MeV with 0.5 x Q 2 and A QC d 
(5 flavors) = 73 MeV with 2.0 x Q 2 for ±er excursions 
around the mean dilepton px values. Since the initial 
and final-state shower algorithms are controlled by the 
same QCD evolution equation [67| . the same variations 
of Aqcd and Q 2 scale are used to estimate the effect 
of final-state radiation. The resulting uncertainty for 
modeling of the initial and final-state radiation is 0.09% 
for CDF and 0.15% for DO. The correction algorithm 
does not distinguish between "soft" (out of cone) and 
"hard" (separate jet) radiation, and there is therefore 
some overlap between the uncertainty on m t for the out- 
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of-cone jet energy correction and for gluon radiation. 
There is also some overlap between the uncertainty for 
initial and final-state radiation and the uncertainty on 
higher-order QCD corrections for high-p^ radiation. 

(iv) Higher-order QCD corrections 

Higher-order QCD corrections to tt production are not 
used for Run I measurements as only LO generators 
were available at that time. DO measures higher- 
order jet-modeling uncertainties in Run II by comparing 
nit extracted with ALPGEN and herwig for evolution 
and fragmentation to the value obtained from events 
generated with MC@NLO [1], which uses herwig parton 
showering with a NLO model for the hard-scattering 
process. This component of uncertainty also includes (for 
DO) the uncertainty from the fraction of quark-antiquark 
to gluon-gluon contributions to the initial state. CDF 
also studies differences in nit using MC@NLO, and finds 
that the uncertainties in distributions in the number 
of jets and the transverse momentum of the tt system 
overlap with the uncertainty from initial and final-state 
radiation. Future measurements of m t are expected to 
treat these uncertainties separately. The uncertainty 
on m t from higher order contributions and initial-state 
qq/gg ratio is 0.14% for DO. 

(v) Jet hadronization and underlying event 

In Run I. CDF measured the uncertainty in the 
model for parton showering and hadronization and the 
underlying event by comparing the value of m t based 
on herwig to that on pythia [§3, and DO compared 
HERWIG results to those from isajet [sH . 

In Run II, CDF estimates these uncertainties by 
comparing m t obtained using pythia with Tune A of the 
underlying event model to results from herwig with a 
tuned implementation of the underlying-event generator 
jimmy. DO estimates these uncertainties by comparing 
identical sets of hard-scatter events from ALPGEN coupled 
to HERWIG instead of to pythia. For the uncertainty on 
nit i this corresponds to 0.40% for CDF and 0.33% for 

do! 

(vi) Color reconnection 

There are up to six final-state quarks in tt events, 
in addition to initial and final-state radiation. When 
hadronization and fragmentation occur, there are color 
interactions among these partons and the color-remnants 
of the proton and antiproton. This process is referred to 
as "color reconnection" . It changes the directions and 
distributions of final-state jets |99l. llOOj . which affects 
the reconstructed value of m t [74j . 

The uncertainty on color reconnection was not 
evaluated for Run I because appropriate MC tools were 
not available at that time. Both collaborations estimate 
this effect in Run II by comparing the value of m t 
extracted from ensembles of tt events generated by 
pythia using the difference between two parton shower 
simulations: (i) angular ordering for jet showers (same as 



used in the nominal mt measurements) using the A-pro 
underlying-event model (Tu ne A but updated using the 
"Professor" tuning tool jlQl|), (") ACR-PRO. ACR-PRO 
is identical to A-pro except that it includes color 
reconnection in the model. The resulting uncertainties 
on mt are 0.32% for CDF and 0.16% for DO. 

5. Multiple interactions model 

Monte-Carlo simulated events are overlaid with 
Poisson-distributcd low-pr events (pythia MC events for 
CDF, "zero-bias" data for DO) to simulate the presence 
of simultaneous additional pp interactions. The mean 
number of overlaid events is chosen at the time of event 
generation, but in data, the number of such interactions 
changes with instantaneous luminosity of the Tevatron. 

CDF measures m t as a function of the number of 
multiple interactions, finding a change of 0.07±0.10 GeV 
per primary vertex. For CDF's measurements, the 
average number of primary vertices in data is 2.20 and for 
simulated events it is 1.85, leading to an uncertainty on 
m t of 0.02%. CDF adds to this in quadrature a term to 
cover the difference in jet energy response as a function 
of the number of multiple interactions of 0.06%, giving a 
total uncertainty of 0.06%. 

DO rcweights the simulated events to make the instan- 
taneous luminosity distribution match that in data. The 
resulting uncertainty on mt is 0.03%. 

6. Background from theory 

(i) Higher-order correction for heavy flavor 

DO corrects the leading-log VF+jets cross section from 
ALPGEN to NLO precision before normalizing this 
background to data. This increases the fraction of Wbb 
and Wcc events in VF+jets by a factor of 1.47 ± 0.50. 
CDF normalizes the ly+heavy-flavor jets background to 
data independent of the other components in VF+jets, 
which has a similar effect. The resulting uncertainties on 
m t are 0.11% for CDF and 0.04% for DO. 

(ii) Factorization scale for W +jets 

The transverse momenta of the jets in VF+jets events 
are sensitive to the factorization and rcnormalization 
scales chosen for the calculations. These two scales are 
set equal to each other, with Q 2 = + J2Pt- To 
determine the uncertainty on m t , the scale is changed 
from (Q/2) 2 to (2 x Q) 2 , the MC events regenerated, 
and the mass remeasured. Changing the scale does not 
affect the fraction of VF+jets in the model but does affect 
the transverse momentum distributions of the jets. The 
uncertainties on m t are 0.02% for CDF and 0.09% for 
DO. 

(in) Normalization to predicted cross sections 

CDF divides the background into seven independent 
parts: W+heavy-flavor jets, W + light-flavor jets, single 
top tqb and tb, Z+jets, dibosons (WW, WZ, and ZZ), 
and multijet contributions. This uncertainty term covers 
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the normalization of the components modeled with MC 
simulated events (not multijets). The small backgrounds 
from single-top, Z+jets, and diboson production are 
normalized to NLO calculations. The uncertainties on 
the cross sections are 10% for tqb, 12% for tb, 14% for 
Z+jets, and 10% for dibosons. The VF+jets background 
is normalized to data before implementation of 6-tagging 
using a fit to the distribution for $t in the event. 
The uncertainty on this normalization cannot be easily 
disentangled from the other sources and so it is kept in 
this category. The combined uncertainty on m t from 
these normalizations is 0.09%. 

DO also normalizes single-top, Z+jets, and diboson 
contributions, in all analysis channels, and Drell- 
Yan in the dilepton channel, to next-to-leading order 
cross sect ions , using values from the mcfm event 
generator [102]. The uncertainties on the cross sections 
take into account the uncertainty on PDF and on the 
choice of factorization and renormalization scales, which 
together propagate through to m t an uncertainty of 
0.04%. 

(iv) Background differential distributions 

For CDF, different methods were used to estimate 
the uncertainty due to the overall background shape. 
In the recent lepton+jets analysis, this uncertainty was 
assessed by dividing randomly the background events 
into subsets, building the background likelihood from one 
of the subsets and reconstruct the mt from the second 
subset. In the next step, the difference in mt obtained 
from the second subset and the nominal mt value is 
evaluated. This contributes an uncertainty of 0.03%. 
CDF also estimates an uncertainty from the limited MC 
statistics used to measure the background. This yields 
an additional 0.03% uncertainty on m t . 

For DO, the pr and 77 distributions of jets in VF+jets 
events do not fully reproduce those in data. An 
uncertainty to cover these deviations is based on the 
difference between the model for background and data 
in the 77 distribution of the third jet in three-jet events. 
The resultant uncertainty on m t is 0.09%. 

7. Background based on data 

(i) Normalization to data 

In the lepton+jets, alljets, $r+jets, and decay- 
length channels, backgrounds from multijet events are 
normalized to data. In the lepton+jets analyses 
at DO, the W+jets background model is combined 
with the contribution from multijet events, and both 
arc normalized simultaneously to data, so that their 
uncertainties in normalization arc anticorrelatcd. In 
dilepton analyses at CDF, the Drell-Yan background is 
normalized to data. For the lepton+jets analyses, CDF 
uncertainty on m t from the normalization of the multijet 
backgrounds to data is 0.03% and DO's uncertainty for 
the normalization of VF+jets and multijets to data is 
0.13%. 



(ii) Trigger modeling 

CDF expects a negligible uncertainty on m t from the 
modeling of the trigger. DO simulates the trigger turn- 
on efficiencies for MC events by applying weights as a 
function of the transverse momentum of each object in 
the trigger. The uncertainty is measured by setting all 
the trigger efficiencies to unity and recalculating the value 
oimt, which shifts m t by 0.03%. 

(iii) b-tagging modeling 

CDF applies the 6- tagging algorithm directly to MC 
events, and finds that any difference between the b- 
tagging behavior in MC and data has a negligible impact 
on the measurement of m t . DO applies the 6-tagging 
algorithm directly to MC events for recent Run II 
measurements. Previously 6-tagging was simulated with 
tag-probability, and in Run I, as DO did not have a silicon 
tracker, nonisolated muons were used to identify 6-jets. 
The tagging efficiency for simulated events is made to 
match that in data by randomly dropping b tags for b 
and c jets, while assigning a per-jet weight for tagging 
light-flavor jets as b jets. The uncertainties for these 
corrections are determined by shifting the efficiencies for 
tagging b and c jets by 5% and by 20% for light jets, 
which introduces an uncertainty on m t of 0.06%. 

(iv) Signal fraction for calibration 

DO measures the impact of the uncertainty in the 
ratio of signal to background events, which affects the 
calibration of m t . Changing the signal fraction within 
uncertainty results in an uncertainty on m t of 0.06%. 

(v) Impact of multijet background on the calibration 

Multijet background events are not used in DO samples 
that determine the calibration of m t for the lepton+jets 
measurement since the background probability for such 
events is much larger than the signal probability. The 
assumption that this has a small effect on mt is tested by 
selecting a multijet-enriched sample of events from data 
(by inverting the lepton isolation criteria) and adding 
these events when deriving the calibration. Applying this 
alternative calibration to data indicates that mt can shift 
by an uncertainty of 0.08%. 

8. Calibration method 

Monte Carlo ti ensembles are generated at different 
values of input m t (m t = 165, 170, 172.5, 175^180 GcV) 
and calibrations relate the input masses for tt events to 
the extracted masses using a straight line. For some 
of the m t measurements, there is an additional in-situ 
calibration of the JES to the light quarks in Vy-boson 
decay, which is then applied to alljets. The uncertainties 
from both calibrations are propagated to the uncertainty 
on m t , which for CDF are 0.04% and 0.05% respectively, 
giving a total of 0.06%. For DO, the uncertainty on m t 
is 0.13%. 
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